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Abstract. The goal of the paper is to develop a Heine-Stieltjes theory for 
univariate linear differential operators of higher order. Namely, for a given 
linear ordinary differential operator = 5Zi=i with polynomial 

coefficients set r = maxi— 1 jj(deg Qi(z) — i). If 0(2:) satisfies the conditions: 
i) r > and ii) deg(3fc(z) = k + r we call it a non-degenerate higher Lame 
operator. Following the classical approach of E. Heine and T. Stieltjes, see I18| . 
[41] we study the multiparameter spectral problem of finding all polynomials 
V{z) of degree at most r such that the equation: 

>}{z)S(z) + V(z)S{z) =0 

has for a given positive integer n a polynomial solution S{z) of degree n. We 
show that under some mild non-degeneracy assumptions there exist exactly 
("n'^) ^^'^^ polynomials Vn^i(z) whose corresponding eigenpolynomials Sn^i{z) 
are of degree n. We generalize a number of well-known results in this area and 
discuss occurring degeneracies. 
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1. Introduction and main results 

The algebraic form of the classical Lame equation, [U], ch. 23, was introduced 
by Gabriel Lame in 1830's in connection with the separation of variables in the 
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Laplace equation in M with respect to elliptic coordinates. It has the form: 

where Qi (z) is a real polynomial of degree I with all real and distinct roots, and V{z) 
is a polynomial of degree at most 1 — 2 whose choice depends on what type of solution 
to (jl.ip we are looking for. In the second half of the 19-th century several celebrated 
mathematicians including M. Bocher, E. Heine, F. Klein, T. Stieltjes studied the 
number and different properties of the so-called Lame polynomials of a given degree 
and certain kind. (They are also called Lame solutions of a certain kind.) Such 
solutions to ()1.1|) exist for certain choices of V{z) and are characterized by the 
property that their logarithmic derivative is a rational function. For a given Q{z) 
of degree I > 2 with simple roots there exist 2' different kinds of Lame polynomials 
depending on whether this solution is smooth at a given root of Q{z) or has there 
a square root singularity, see details in [32j and [17j . (An excellent modern study 
of these questions can be found in [20].) In what follows we will concentrate on the 
usual polynomial solutions of (jl.ip and its various modifications. 

A generalized Lame equation, see |47| is the second order differential equation 
given by 

Q2{z)^ + Qi{z)^ + V{z)S = Q, (1.2) 
dz 

where Q2iz) is a complex polynomial of degree I and Qi{z) is a complex polynom of 
degree at most I — 1. The special case I = 3 is widely known as the Heun equation. 

The next fundamental proposition announced in [18] and provided there with 
not a quite satisfactory proof was undoubtedly the starting point of the classical 
Heine-Stieltjes theory. 

Theorem 1 (Heine). If the coefficients of Q2{z) and Qi{z) are algebraically in- 
dependent, i.e. they do not satisfy an algebraic equation with integer coefficients 
then for any integer n > there exists exactly ("^^~^) polynomials V{z) of degree 
exactly (/ — 2) such that the equation has and unique (up to a constant factor) 

polynomial solution S of degree exactly n. 

Remark 1. Notice that throughout this paper we count polynomials V{z) individ- 
ually and polynomials S{z) projectively, i.e. up to a constant factor. 

Later on a physically important and directly related to the original (jl.ip special 
case of (|1.2p when Q2{z) and Qi{z) have all real, simple and interlacing zeros and 
the same sign of the leading coefficients was considered separately by T. Stieltjes 
and his followers. The equation can be then written as follows: 

j=i 

with ai < a2 < ■ ■ ■ < Oil real and (3i,...,Pi positive. In particular, the next 
proposition was proved. 

Theorem 2 (Stieltjes- Van Vleck-Bocher [H], [IS], [5] and [12]). Under the as- 
sumptions of U.3]) and for any integer n > 

(1) there exist exactly ("^^~^) distinct polynomials V of degree (1 — 2) such that 
the equation has a polynomial solution S of degree exactly n. 

(2) each root of each V and S is real and simple, and belongs to the interval 
{ai,ai). 

(3) none of the roots of S can coincide with some of ai 's. Moreover, ("^^^^) 
polynomials S are in 1-1-correspondence with ("~^^~^) possible ways to dis- 
tribute n points into the (^ — 1) open intervals (01,02), (02, c^s),- • • ; (cti-ijCti) 
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The polynomials V and the corresponding polynomial solutions S of the equa- 
tion p.2p (or, equivalently, of (|1.3p ) are called Van Vleck and Stieltjes (or Heine- 
Stieltjes) polynomials resp. 

The case when ai's and/or /3j's are complex is substantially less studied, see [SS] 
and [U]. One nice result in this set-up is as follows, see [HI]- 

Theorem 3 (Polya). // in the notation of il.3\} all ai 's are complex and all fij 's 
are positive that all the roots of each V and S belong to the convex hull Convg^ of 
the set of roots (ai, . . . , a;) of Q2{z). 

Remark 2. The situation when all the residues Pj are negative (for example, 
Qi{z) = —Q2{z)) or have different signs seems to differ drastically from the latter 
case, see e.g. ^44^ and [H]. Further interesting results on the distribution of the 
zeros of Van Vleck and Stieltjes polynomials under weaker assumptions on a^'s and 
/3/s were obtained in [2T]. p]. [2]. [48|. 



In the present article we extend the above three fundamental results on general- 
ized Lame equations of the second order to the case of higher orders and/or complex 
coefficients. Namely, consider an arbitrary linear ordinary differential operator 



with polynomial coefficients. The number r = TOaa;i=i_..._fc(deg (5i(z) — i) will be 
called the Fuchs index of 0(2:). The operator 0(0) is called a higher Lame operator 
if its Fuchs index r is non-negative. In the case of the vanishing Fuchs index i)(z) is 
usually called exactly solvable in the physics literature, see [13] ■ This case is also of 
the special interest in connection with the classical Bochner-Krall problem in the 
theory of orthogonal polynomials. 

The operator i)(z) is called non-degenerate if degQk{z) = fc + r. Notice, that 
non-degeneracy of 5(z) is a quite natural condition equivalent to the requirement 
that D(z) has either a regular or a regular singular point at 00. 

Given a higher Lame operator 5(z) consider the multiparameter spectral problem 
as follows. 

Problem. For each positive integer n find all polynomials V{z) of degree at most 
r such that the equation 



has a polynomial solution S{z) of degree n. 

Following the classical terminology we call (II. 5|) a higher Heine- Stieltjes spec- 
tral problem, ^i^) is called a higher Van Vleck polynomial, and the corresponding 
polynomial S{z) is called a higher Stieltjes polynomial. Below we will often skip 
mentioning 'higher'. 

Remark 3. Obviously, any differential operator (|1.4p has either a non-negative or 
a negative Fuchs index. In the latter case it can be easily transformed into the 
operator with a non- negative Fuchs index by the change of variable V = \- Notice 
also that the condition of non-degeneracy is generically satisfied. In what follows 
we will always assume wlog that the leading coefficient of such an operator is a 
monic polynomial. 




k 



(1.4) 



i){z)S{z) + V{z)S{z)=Q 



(1.5) 
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1.1. Generalizations of Heine's theorem, degeneracies and nonresonance 
condition. Wc start with a number of generalizations of Heine's theorem [T] Fol- 
lowing Heine's original proof one can obtain the following straightforward general- 
ization. 

Theorem 4. For any non-degenerate higher Lame operator X){z) with algebraically 
independent coefficients of its polynomial coefficients Qi{z), i = l,...,k and for 
any n > there exist exactly ("^'") distinct Van Vleck polynomials V{z) 's whose 
corresponding Stieltjes polynomials S{z) 's are unique (up to a constant factor) and 
of degree n. 

Our next result obtained by a linear- algebraic interpretation of ()1.5p has no 
genericity assumptions and is crucial in the problem of existence of solutions of 
(|L5l) . comp. [271, Problem 1. 

Theorem 5. For any non-degenerate operator with a Fuchs index r > and 
any positive integer n the total number of Van Vleck polynomials V{z) (counted 
with natural multiplicities) having a Stieltjes polynomial S{z) of degree less than or 
equal to n equals ("^_]^^^)- 

Remark 4. Note that in Theorem O we do not require that there is a unique 
(up to constant factor) Stieltjes polynomial corresponding to a given Van Vleck 
polynomial. For the exact description of the notion of the natural multiplicity of a 
Van Vleck polynomial which is rather lengthy consult Definition [5] in § [21 

On degeneracies. Notice that Theorem [J claims that a generic operator c)(z) 
has for any positive n exactly ("^'^) distinct Van Vleck polynomials each of which 
has a unique Stieltjes polynomial and this polynomial is of degree exactly n. The 
question about possible degeneracies occurring in Problems ()1.2p and ()1.5p if we 
drop the genericity assumptions on is quite delicate. Not only Van Vleck 
polynomials can attain a nontrivial multiplicity as well as more than 1-dimensional 
linear space of Stieltjes polynomials but there are examples when there are no 
Stieltjes polynomials of some degree. In particular, for any polynomial Q{z) of 
degree I no choice of a polynomial V{z) of degree at most I — 2 will supply the 
equation 

Q{z)^-Q'iz)^ + V{z)S = 
dz'^ dz 

with a polynomial solution S of degree ^ -f- 1. (This follows from the Proposition 
5 and Lemma 4 of ^14].) The fact that (jl.2p can admit families (linear spaces of 
dimension at least 2) of polynomial solutions S corresponding to one and the same 
V was already mentioned by Heine in his original proof. More exact information is 
available nowadays. For example, a result of Varchenko-Scherbak gives necessary 
and sufficient condition for a Fuchsian second order equation to have 2 independent 
polynomial solutions, see [21] and [2]. Finally, high multiplicity of Van Vleck 
polynomials occur, for example, in the case Q2{z) = z', Qi{z) = 0. Then one can 
easily show that for all n > 2 there exists just one and only polynomial Vn{z) = 
—n{n — l)z'~^ solving the above problem; its corresponding Stieltjes polynomial 
equals Sn{z) = z^. The multiplicity of the latter Van Vleck polynomial Ki(z) is 

To formulate necessary and sufficient conditions under which the conclusion of 
Heine's theorem holds for all positive integers n is apparently an impossible task. 
Heine himself mentions that for the validity of his result for a given fixed positive 
integer n one has to avoid a certain discriminantal hypersurface (similar to the 
usual discriminant of univariate polynomials) which is given by an equation with 
integer coefficients but this equation is difficult to obtain explicitly. 
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Below we formulate a simple sufficient condition which allows us to avoid many 
of the above degeneracies and guarantees the existence of Stieltjes polynomials 
of a given degree. Namely, consider an arbitrary non-degenerate operator £)(z) 
of the form ()1.4|) with the Fuchs index r. Denote by Ak, Ak-i, Ai the coeffi- 
cients at the highest possible degrees k + r, k + r — 1, ...,r + 1 in the polynomials 
Qk{z),Qk^i{z), ...,Qi(z) resp. (Notice that any subset of Aj^s can vanish but 
Ak ^ due to the non-degeneracy of f (z)-) In what follows we will often use the 
notation 

(i). =j(j-i)(i-2)...(j-i + i), 

where j is a non-negative and i is a positive integer. In case j = i one has — jl 
and in case j < i one gets = 0. For any non-negative n we call by the n-th 
diagonal coefficient L„ the expression: 

L„ = in)kAk + {n)k-iAk-i + .... + (n)iAi. (1.6) 

Proposition 1. // in the above notation and for a given positive integer n the n-th 
nonresonance condition 

L„^L,, j = 0,l,-,n-l (1.7) 

holds then there exist Van Vleck polynomials which possess Stieltjes polynomials of 
degree exactly n and no other Stieltjes polynomials of degree smaller than n. In 
this case the total number of such Van Vleck polynomials (counted with natural 
multiplicities) equals ("^'^) • 

Remark 5. The above nonresonance condition is quite natural. It says that if the 
equation p.Sp has a polynomial solution of degree n then it has no polynomial 
solutions of smaller degrees. Another way to express this fact is that if the indicial 
equation of (|1.5p at oo has — n as its root then it has no roots among non-positive 
integers 0, —1, —2, 1 — n, see e.g. [32], ch. V. 

Explicit formula (|1.6p for L„ immediately shows that Theorem [5] and Proposi- 
tion [1] are valid for any non-degenerate d{z) and all sufficiently large n. 

Corollary 1. For any non- degenerate higher Lame operator d{z) and all sufficiently 
large n the n-th nonresonance condition holds. In particular, for any problem lll.5\) 
there exist and finitely many (up to a scalar multiple) Stieltjes polynomials of any 
sufficiently large degree. 

Remark 6. Notice that for an arbitrary non-degenerate operator d{z) and a given 
integer n it is difficult to find explicitly all Van Vleck polynomials which possess 
a Stieltjes polynomial of degree at most n. By this we mean that in order to do 
this one has, in general, to solve an overdetermined system of algebraic equations 
in the coefficients of V since the set of Van Vleck polynomials under consideration 
is not a complete intersection. (This system of determinantal equations contains 
many more equations than variables.) However, one consequence of Heine's way to 
prove his Theorem [T] is as follows. For a given non-degenerate operator d{z) with 
Fuchs index r and a positive integer n denote by 2J„ C Polr the set of all its Van 
Vleck polynomials possessing a Stieltjes polynomial of degree exactly n. 

Theorem 6. // in the above notation the n-th nonresonance condition fi. 7| ) holds 
then QJ„ is a complete intersection and the corresponding system of equations can 
be given explicitly in each specific case. 

Explicit example of the defining system of r algebraic equations in r variables 
can be found in § [21 see Example [Tj In purely linear algebraic setting this result 
and further information about relevant discriminants can be found in [5U] . 
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1.2. Generalizations of Stieltjes's theorem. We continue with a conceptually 
new generalization of Theorem [21 It was conjectured by the present author after 
extensive computer experiments and was later proved by P. Branden. In the present 
paper we only announce his result and its corollaries since it requires a large amount 
of additional information and techniques. The actual proof will be published by its 
author elsewhere. 

Definition 1. A differential operator d{z) — X]i=m Qii^)'IF' 1 < < ^ where all 
(5i(z)'s are polynomials with real coefficients is called a strict hyperbolicity preserver 
if for any real polynomial P{z) with all real and simple roots the image d{P{z)) 
either vanishes identically or is a polynomial with only real and simple roots. 

Theorem 7. For any strict hyperbolicity preserving non-degenerate Lame operator 
t)(z) with the Fuchs index r as above and any integer n > m 

(1) there exist exactly {"^^) distinct polynomials V{z) of degree exactly r such 
that the equation U.5\) has a polynomial solution S{z) of degree exactly n. 

(2) all roots of each such V{z) and S{z) are real, simple, coprime. 

(3) {"^^) polynomials S{z) are in 1- 1-correspondence with ("^'^) possible ar- 
rangements of r real roots of polynomials V(z) and n real roots of the cor- 
responding polynomials S{z). 

Using Theorem [5] one immediately sees that the latter result describes the set of 
all possible pairs (V, S) with m < n — degS for any hyperbolicity preserver d{z). 

Remark 7. The interested reader can check that the sum of the first two terms in 
(|1.3p is indeed a strict hyperbolicity preserver. It looks very tempting and important 
to find an analog of the electrostatic interpretation of the roots of classical Heine- 
Stieltjes and classical Van Vleck polynomials (alias 'Bethe ansatz') in the case of 
higher Heine-Stieltjes and Van Vleck polynomials, comp. [28] . 

Remark 8. Notice that the converse to the above theorem is false. Namely, one can 
show that the exactly solvable operator D{z){f) = f + z{z + l)f" has all hyperbolic 
eigenpolynomials but is not a hyperbolicity preserver. 

A straight-forward application of Theorem [7] to differential operators of order 2 
gives the following. Consider a differential equation 

ri(- - IK- - + = 0' ^'-^^ 

i=l j=l j^i 

where 2 < k < I, ai < a2 < ■ ■ ■ < ai and . . . , /3; are positive. 
Corollary 2. Under the assumptions of m.8\) and for any n > k ~ I 

(1) there exist exactly polynomials V{z) of degree {I — k) such that the 
equation hl.^) has a polynomial solution S{z) of degree exactly n. 

(2) all roots of each V{z) and S{z) are real, simple, coprime and belong to the 
interval (ai, ai). 

(3) ("^^~'^) polynomials S{z) are in 1- 1-correspondence with possible 
arrangements of {I — k) real roots of a polynomial V{z) and n real roots of 
the corresponding polynomial S{z) on the interval {ai,ai). 

It seems that Theorem[7]and Corollary [5] give a new interpretation of Theorem[5] 
even in the classical case (|1.3p . However the following two statements proven by 
G. Shah explain this mystery, see Theorem 3 of [35] and Theorem 3 of [38j . 
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Proposition 2. Under the assumptions of Theorem\^the roots of any Van Vleck 
polynomial V{z) and its corresponding S{z) are coprime. 

Moreover, 

Proposition 3. Ifvi < V2 < ■ ■ ■ < f-r, f ~ 1 — 2 denote the roots of some Van Vleck 
polynomial V{z) in the classical situation then for each i — 2. . . . , Z — 1 the 

interval (wijOfi+i) contains no roots of the corresponding S{z). Therefore, for each 
polynomial S{z) the distribution of its n roots into 1) intervals (ai, 02), ■ ■ ■ , a/) 
coincides with the distribution of these roots defined by the roots of its Van Vleck 
polynomial V{z). 

Remark 9. Note that we do not claim Vi < Ofi+i, i.e. the endpoints of the interval 
{vi,ai+i) can be placed in the wrong order or can coincide. 

1.3. Generalizations of Polya's theorem. 

We start with a simple-minded statement of Polya's theorem [31 [31] . 

Theorem 8. If the zeros ai, . . . , a/ in il.8\} are complex and the constants . . . , /?/ 
are non-negative then all the roots of V 's and S 's lie in the (closed) convex hull 
ConvQf, of the roots (ai, . . . ,ai) of the polynomial Qk{z). 

The next Theorem[H|is far more general. It shows that Theorem[H]is asymptoti- 
cally true for any non-degenerate Lame operator. The question for which operators 
X){z) the roots of all its Van Vleck and Stieltjes polynomials lie exactly (and not 
just asymptotically) in the convex hull of its leading coefhcient seems to be very 
hard even in the classical case of the equation (|1.2p . 

Theorem 9. For any non- degenerate higher Lame operator X){z) and any e > 
there exists a positive integer such that the zeros of all Van Vleck polynomials 
V{z) possessing a Heine- Stieltjes polynomial S{z) of degree n > and well as all 
zeros of these Stieltjes polynomials belong to Convq^ . Here Convg^, is the convex 
hull of all zeros of the leading coefficient Qk and ConVg^ is its e-neighborhood in 
the usual Euclidean distance on C. 

The latter theorem is closely related to the next somewhat simpler localization 
result having independent interest. 

Proposition 4. For any non- degenerate higher Lame operator d{z) there exist a 
positive integer Nq and a positive number Rq such that all zeros of all Van Vleck 
polynomials V{z) possessing a Stieltjes polynomial S{z) of degree n > Nq as well 
as all zeros of these Stieltjes polynomials lie in the disk \z\ < Rq. 

Remark 10. Notice that the roots of absolutely all Van Vleck polynomials (and 
not just those whose Stieltjes polynomials are of sufficiently large degree) of any 
non-degenerate higher Lame operator 0(2:) lie in some disk. But this is no longer 
true for Stieltjes polynomials. If the set of all Stieltjes polynomials is discrete (up 
to a scalar multiple) then their roots are bounded. But as soon as some Van Vleck 
polynomial admits an at least 2-dimensional linear space of Stieltjes polynomials 
then these roots become unbounded for obvious reasons. However for sufficiently 
large n no Van Vleck polynomial admits such families, see Corollary [T] and the 
localization result holds. 

Remark 11. Similar and stronger localization results with explicit constants and 
degree bounds were independently obtained by J. Borcea (private communication). 

Let us now show a typical behavior of the zeros of Van Vleck polynomials and the 
corresponding Stieltjes polynomials obtained in numerical experiments. Below we 
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consider as an example the operator = Q{z)^ with Q{z) — [z^ -\- l){z — il — 
2){z + 2/ - 3). For n = 24 we calculate ah 25 pairs {V, S) with deg S = 24. (Notice 
that V in this case is linear.) The asymptotic behavior of the union of zeros of all 
Van Vleck polynomials whose Stieltjes polynomials have a given degree n when n — > 
oo as well as the asymptotics of the zeros of subsequences of Stieltjes polynomials 
of increasing degrees whose corresponding (monic) Van Vleck polynomials have a 
limit seems to be an extremely rich and interesting topic, see first steps in [40j. 



1 • 

-1 • • • 

-2 • 

-3L . . ^ 



Figure 1. Zeros of 25 different linear Van Vleck polynomials 
whose Stieltjes polynomials are of degree 24. Four average size 
dots are the zeros of Q{z) = {z^ + l)(z - 3/ - 2){z + 21-3) and 
large dots are the zeros of different V{z). 



Figure 2. Zeros of 25 different Stieltjes polynomials of degree 24 
for the above 5(2). The small dots on each picture are the 24 zeros 
of S{z); 4 average size dots are the zeros of Q{z) and the single 
large dot is the (only) zero of the corresponding V{z). 



Some literature. Let us mention a few relatively recent references on (general- 
ized) Lame equation. Being an object of substantial physical and mathematical 
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importance it, in particular, gives an example of an equation whose monodromy 
group can be analyzed in details, see [7], [Il]-[2Dj. It is also closely related to the 
so-called quasi-exact solvability and integrable models, see [H]. Theory of multipa- 
rameter spectral problems originating from Heine-Stieltjes pioneering studies was 
developped in sixtees, see e.g. [46 and references therein. Recently the interest 
to Heine-Stieltjes polynomials has stimulated by an unexpected extension of the 
Bethe ansatz in representation theory, see [331 1 i [31] , [31] and a further series of 
article by A. Varchenko and his coauthors, e.g. [30]. Starting with [35] a substantial 
progress has been made in the understanding of the asymptotics of the root distri- 
butions for these polynomials when either I — > oo (thermodynamic asymptotics) or 
n — > oo (semi-classical asymptotics), see |10| . and [S]. Asymptotic root distribution 
for the eigenpolynomials of non-degenerate exactly solvable operators was studied 
in [24j and [H]. Interesting preliminary results of a similar flavor in the case of 
degenerate exactly solvable operators were very recently obtained by T. Bergkvist, 

Acknoledgements. I am very grateful to my former collaborator and colleague 
G. Masson for the pioneering numerical experiments in 1999. His wild guesses 
gave birth to a vast project on asymptotics for polynomial solutions to linear or- 
dinary differential equations depending on parameter(s) which occupies me since 
then. I want to thank P. Branden for sharing his proof of Theorem [7| with me in 
a private communication. Sincere thanks go to R. B0gvad, J. Borcea, I. Scherbak, 
A. Varchenko and, especially, to A. Martinez-Finkelshtein for many useful discus- 
sions of the area and their interest in my work. Finally, I owe a great deal to the 
Wolfram corporation whose package Mathematica although quite expensive and 
not quite reliable was indispensable in doing the actual fiaOrjiianKa. 



2. Proof of generalized Heine's theorems 

We start with Theorem (see Introduction) . For this we need a detailed de- 
scription of the action of a non-degenerate operator d{z) on the linear space Poln 
of all univariate polynomials of degree at most n. 



Proof. Substituting V{z) = VrZ'^ +Vr-iz^~^ + ■ ■ .+wo and S{z) = SnZ^ + Sn-iz'^~'^ + 
... -I- So in (|1.5p we get the following system of [n + r-\-l) equations of a band shape 
(i.e. only a fixed and independent of n number of diagonals is non- vanishing in this 
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system): 

= Sn{Vr + Ln^n+r)', 

= S„(Wr-l + Ln^n+r-l) + Sn-l{Vr + -^n- l,n+r-l ) ; 

= Sn{Vr-2 + l{Vr-l + Ln-ljn+r-2) + Sn-2{Vr + Ln-2,n+r-2)] 

= S„(wo + Ln,n) + Sn-l(wi + Ln-l,n) + ••• + S„- 

= S„L„^„_i + S„_i(wo + -^n-l,n-l) + ■■■ + S„-r-l('yr + in-r- l,n- 1 ) ; 
< = SnLn,n-2 + + S2(Wr + -^2,r) + ••• + S„-r-2(''^r + iri-'r-2.n-2) ; 

= SnLn.r + Sn-l-^n-l,r + S„_2in-2,r + ... + So(Ur + Lq^t)', 

= SnLn,r-l + Sn-l-^n-l,r-l + Sn-2-^n-2,r-l + ■•. + So(''^r-l + -^0,l); 

= s„L„4 + s„-iL„_i,i + ... + .si(i;o + Li.i) + so(«i + -^0,1); 
^0 = s„L„.o + s„-iL„_i.o + ... + siLo,i + so(wo + -^o.o). 

(2.1) 

Here Lp ^ is a polynomial which expresses the coefficient containing Sp at the 
power z"^ in X]i=i Qi{z)S^^K Obviously, it is linear in the coefficients ofQk{z), Qi{z) 
and is explicitly given by the relation 

k 

^p,q — ^ ^ (p)r^r,<j— p+r; 
r=l 

where Aj.,g_p+r is the coefficient at 2'"^'+'' in Qr{z). In the notation used in the 
definition p.6p we have Lm,m+r = ^m,rn = 0, We use the convention that 

Lp^q vanishes outside the admissible range of indices and, therefore many of the 
above coefficients Lp^q are in fact equal to 0. (In the system (|2.ip we assumed 
that n > r for simplicity.) Notice that all equations in (|2.ip depend linearly on 
the variables Vr,...,VQ and s„,...,so as well as on the coefficients of polynomials 
Qi{z), i — 1, . . . ,/c. Note additionally, that (|2.ip is lower-triangular w.r.t the co- 
efficients Sn, So which allows us to perform the following important elimination. 
Let us enumerate the equations of (|2.ip from to n -I- r assigning the number j to 
the equation describing the vanishing of the coefficient at the power . Then 

if Lj, = Ln_n+r One has that the 0-th equation has a solution s„ = 1 and Vr = 
—Ln^n+r 7^ 0. The next n equations are triangular w.r.t the coefficients s„, sq, i.e. 
j-th equation in this group contains only the variables s„, s„_i, . . . , s„_j (among 
all Sj's) along with other types of variables. Thus under the assumption that all 
the diagonal terms Vr + Ln-i^n+r-i = ^n-i — Ln, i = 0, 1, ...,n are nonvanishing 
we can express all Sn~i, i ~ 0,1,..., n consecutively as rational functions of the 
remaining variables and get the reduced system of r rational equations containing 
only (fr-i, . . . , vq) as unknowns. Notice that in view of u,. = —Ln.n+r the non- 
vanishing of the diagonal entries Vr + Ln-i,n+r-ii * = 0, 1, coincides exactly 
with the nonresonance condition (|1.7p . 

Cleaning the common denominators we get a reduced system of polynomial equa- 
tions. We show now that this polynomial system is quasi-homogeneous in the vari- 
ables Vj with the quasi-homogeneous weights w{vj) given by w{vj) = r — j . Thus 
using the weighted-homogeneous version of the Bezout theorem, see e.g [13] we 
get that if the system under consideration defines a complete intersection, i.e. has 
only isolated solutions then their number (counted with multiplicities) equals ("^'') 
in the corresponding weighted projective space. To check the quasi-homogenuity 
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note that the standard action of C* on the set of roots of the polynomial V{z) by 
simultaneous multiplication assigns the weight r ~ j to its coefficient Vj. These 
weights are still valid in the reduced system with the variables s„, sq eliminated. 
Finally, we have to show that if the coefficients of Qk{z), Qi{z) are algebraically 
independent then the eliminated system has exactly ("^'^) simple solutions. Indeed, 
consider the linear space EQ of all systems of r quasi-homogeneous equations in 
the variables (vr, ■■■,vo) with the weights w{vj) — r — j and where the z-th equation 
is weighted- homogeneous of degree n + i. We equip this space with the standard 
monomial basis. 

To accomplish the proof of Theorem |4] we need two additional standard facts. 

Lemma 1. The discriminant Discr C EQ (i.e. the set of the coefficients of mono- 
mials in the equations for which the system has at least one solution of multiplicity 
greater than I) is given by an algebraic equation with rational coefficients in the 
standard monomial basis of EQ. 

Proof See [H], ch. 13. □ 

Consider some linear parameter space A with a chosen basis. Assume there is a 
rational map: $ : A — > EQ where each coordinate in the standard monomial basis 
of EQ is given by a rational function with rational coefficients w.r.t to the chosen 
basis in A. The next statement is obvious. 

Lemma 2. In the above notation the pullback of {Discr) in A either a) coin- 
cides with the whole K or b) is given in the chosen basis by an algebraic equation 
with rational coefficients. 



It remains to show that there are some values of the coefficients of the polyno- 
mials Qk{z), Qi{z) for which there are exactly ("^'^) distinct solutions of (|1.5p . 
Here we are not able to follow the nice inductive argument of [TB], see also the 
last paragraph in § |31 Heine's proof does not generalize immediately to higher 
order equations. Instead we can, for example, invoke Theorem [7] whose proof is 
completely independent of the present arguments. It claims, in particular, that 
for any strict hyperbolicity preserver of the form X){z) — X]i=m ^'^'i ^"^y 

n > m there exist exactly ("^'^) pairs {V, S) . One can additionally choose such a 
hyperbolicity preserver with to = 1 and therefore get a necessary example of an 
operator with given k and r such that for any n > 1 it has exactly the maximal 
number of pairs (F, S) . □ 

To settle Theorem [S] (see Introduction) let us first reinterpret Problem p.5p in 
linear algebraic terms. 

2.1. On eigenvalues for rectangular matrices. We start with the following 
natural question. 

Problem. Given a l)-tuple of (toi x TO2)-matrices A, i?; where to i < m2 

describe the set of all values of parameters Ai, ...A; for which the rank of the linear 
combination A + XiBi + ... -\- XiBi is less than toi i.e. when the linear system 
V * {A -\- XiBi + ... + XiBi) = has a nontrivial (left) solution v ^ which we call 
an eigenvector of A wrt the linear span of Bi, ...,Bi. 

Let A4mi.m2 denote the linear space of all (toi x TO2)-matrices with complex 
entries. Below we will consider l-tuples of (mi x m2)-matrices Bi, Bi which are 
linearly independent in AArni,m2 a-nd denote their linear span hy C = C{Bi, ...,Bi). 
Given a matrix pencil V ~ A + C where A 6 Mmi,m2 denote hy £-p C V its 
eigenvalue locus, i.e. the set of matrices in V whose rank is less than the maximal 
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one. Denote by A^^ C Mmi,m2 the set of all (mi x 7712) matrices with positive 
corank, i.e whose rank is less than mi. Its co-dimension equals m2 — mi + 1 and 
its degree as an algebraic variety equals („™ii), see [12], Prop. 2.15. Consider 
the natural left-right action of the group GLmi x GLm2 on A4mi,m2j where GLm^ 
(resp. GLm.-^) acts on (mi x TO2)-matrices by the left (resp. right) multiplication. 
This action on M-rni,m2 has finitely many orbits, each orbit being the set of all 
matrices of a given (co)rank, see e.g. |4;, ch.l §2. Notice that due to the well- 
known formula of the product of coranks the codimension of the set of matrices 
of rank < r equals (toi — r){m2 — r). Obviously, for any pencil V one has that 
the eigenvalue locus coincides with = A^^ n "P. Thus for a generic pencil V of 
dimension I the eigenvalue locus £-p is a subvariety of V of codimension m2 — mi -I- 1 
if I > 7712 — "nil + 1 and it is empty otherwise. The most interesting situation for 
applications occurs when I — m2 — mi -I- 1 in which case E-p is generically a finite set. 
From now on let us assume that I — m2 — mi -f 1. Denoting as above by C the linear 
span of Bi,...,Bi we say that C is transversal to if the intersection C D 
is finite and non-transversal to otherwise. Notice that due to homogeneity of 
Ai^ any (m2 — mi -|- l)-dimensional linear subspace C transversal to it intersects 
only at and that the multiplicity of this intersection at equals (^™ii). 
We start with the following obvious statement which will later imply Theorem[5l 

Lemma 3. // (m2 — mi -t- 1)- dimensional linear space C is tranversal to then 
for any matrix A G A4mi,m2 the eigenvalue locus £-p of the pencil V = A + C 
consists of exactly („™^i) points counted with multiplicitites. 

Remark 12. Notice that since A^^ C Al(mi,m2) is an incomplete intersection then 
in order to explicitly determine the eigenvalue locus of a given matrix A w.r.t. some 
(m2 — mi -|- l)-dimensional linear subspace C C A^(mi,m2) one has to solve an 
overdetermined system of equations describing the vanishing of all maximal 
minors of a (mi x m2)-matrix depending on parameters. Fortunately, in our main 
application, i.e. for the multi-parameter spectral problem (jl.Sp we encounter only 
'triangular' rectangular matrices (i.e. with the left-lower corner vanishing) for which 
the determination of the eigenvalue locus often reduces to a complete intersection, 
see proof of Theorem 2] and Theorem [S] 

Let us explain how Lemma[3]implies Theorem^ Namely, given a non-degenerate 
operator X){z) in order to find all its Van Vleck polynomials having (at least one) 
Stieltjes polynomial of degree at most n we need to study the action of 0(2) on the 
linear space Poln of all univariate polynomials of degree at most n. If X){z) has the 
Fuchs index r then 5(2;) maps Poln to Poln+r- Using the standard monomial basis 
1, z, z^, in Poll we get that if n > fc = ord(5(z)) then the action of X){z) in this 
basis is represented by a 'triangular' band (n-l-1) x (ri-l-r-l-l)-matrix >lo(z)^„ with at 
most r + fc non-vanishing diagonals. Here 'triangular' means that all entries j of 
^o(z),n with i < j vanish. Denote by Is, s = 0, the (n-l-1) x (n -I- r -I- l)-matrix 
whose entries are given by aij = if i — j 7^ s and 1 otherwise. Denote by £ the 
linear span of /q, Ir and notice that £ is transversal to J\A^ C Af„+i,„+r+i since 
any matrix belonging to the pencil £ and different from has full rank. 

Notice that adding an arbitrary polynomial V{x) = VrZ^ + Vr-iz^^^ + ... + vq 
of degree at most r to 5(z) corresponds on the matrix level to adding of the linear 
combination Vrlo + Wr-i^i + ... + v^Ir to the initial matrix ^^(z) „. The existence 
of a non-trivial Stieltjes polynomial of degree at most n corresponds to the fact 
that the matrix Agj^) „ -I- Vrlo + Vr-ili + ... + vqIt has a non-trivial (left) kernel. 
Thus, for a given non-degenerate operator d{z) the problem of finding all Van 
Vleck polynomials whose Stieltjes polynomials are of degree at most n is exactly 
equivalent to the determination of all the eigenvalues of its matrix Agj^) „ w.r.t. 
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the Hnear space £ in the above-mentioned sense. Lemma [3] has a simple analog for 
'triangular' rectangular matrices which is equivalent to Theorem [5l 

Namely, denote by TM{mi,m2) C M{mi,m2), mi < m2 the set of all 'tri- 
angular' mi X m2-matrices, i.e. with Uij — for i < j. Let £ C TM{mi, 7112) be the 
linear subspace spanned by all (7712— mi -1-1) possible unit matrices /i, /ma-mi+i £ 
TM{mi,m2). Finally, denote by TM^ C TM{mi,m2) the set of aU 'triangular' 
matrices with positive corank. 

Lemma 4. For any matrix A G TM{nii, m2) the eigenvalue locus S-p of the pencil 
V = A + £, consists of exactly („™li) points counted with multiplicitites. 

Proof. The same as above. □ 



The latter Lemma settles Theorem [51 Let us now prove Proposition [T] (see 
Introduction). 

Proof. In the above notation consider the pencil An(vr,Vr-i, ...,vq) = ^o(z),n + 
Vrlo -\- Vr-ih + ... -|- Voir of {n + 1) X {n + r + l)-matrices. One has 

/L„ + Vr * * * * • • A 

Ln-l + Vr * * * • • 

IL„_2 + Vr * * ■ ■ 

L„_3 +Vr * ■■ 



A„{Vr,Vr-l, ...,Wo) = 



V : : : : :/ 

where * stands for possibly non-vanishing entries. The obvious necessary condition 
for such a matrix to have a positive corank it that one of the elements on the shown 
above main diagonal vanishes, i.e. there exists i = 0, ...,n such that ILi + Vr = 
or, equivalently, Vr = — L.^. Set Vr = — L„ thus 'killing' the entry in the left-upper 
corner. Recall that the n-th nonresonance condition requires that L„ ^jij = 
0, ...,n — 1. Therefore the subtraction of L„ along the main diagonal will keep 
all other diagonal entries except for the left-upper corner non-vanishing. The r- 
dimensional pencil y^„(— L„, w^-i, vq) = ^t)(z).n — L„/o + f r-i-^i + ••■ + voir has 
its 1-st column of the above matrix presentation vanishing and its main diagonal 
being the same for all possible values of tj^-i, wq- Therefore, by Lemma [3] there 
exist exactly {"^^) eigenvalues w^-i, Wr-2j "I'd counted with multiplicities such 
that L„, i;j.-i, Wo) has a positive corank. Finally, notice that since for any 

matrix from the pencil yl„(— L„, t^^-i, wo) its entries along the main diagonal 
except for the left-upper corner are non- vanishing its corank can be at most 1. 
Moreover, when the corank of such a matrix is 1 then the occurring non-trivial 
linear combination of the rows which vanishes must necessarily include the first row 
since the second, the third etc rows are linearly independent for the above reason. 
The coefficients of this linear dependence of rows are exactly the coefficients of 
the corresponding Stieltjes polynomial. The fact that the first row must be in the 
linear dependence means in this language that the leading coefficient of this Stieltjes 
polynomial (which by definition is of degree at most n) must be non- vanishing, i.e. 
this Stieltjes polynomial is of degree exactly n. Proposition [1] is settled. □ 

To prove Theorem [S] (see Introduction) we need to take a more careful look at 
the proof of Theorem U] Namely, consider again the system (|2.ip determining the 
set of all pairs (V, S) where F is a Van Vleck polynomial and S is the corresponding 
Stieltjes polynomial of degree at most n. As in the proof of Theorem |4] we solve 
the 0-th equation in (j2.ip by taking s„ = 1 and Vr = —Ln^n+r ~ ~'^n- (By 
Proposition[T]finding a solution of l|2.ip with s„ = 1 and tv = —Ln.n+r — — L„ leads 
to a pair (F, S) such that V is of degree exactly equal to r and S is of degree exactly 
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equal to n.) Then we express consecutively the variables s„_i, s„_2, sq from the 
next n equations of (|2.ip . The crucial circumstance here is that while doing this we 
only divide by the differences of the form L„ — Lj , j = n — 1, rt — 2, which are 
non-vanishing due to the validity of our nonresonancc condition. Substituting the 
obtained expressions for Sj, j = 0, n in the remaining r equations in (|2.ip we get 
the required eliminated system of algebraic equations on the variables Vr^i, ...,wo 
which proves the required result. □ 
To illustrate the above procedure let us consider a concrete example. 

Example 1. Consider the action of some operator d{z) with the Fuchs index r = 2 
on the space Poh. Its maps Poli to P0I3 and, say, is represented in the monomial 
bases of Poll and P0I3 by the matrix 



(Here we used the notation from the proof of Theorem Since r — 2 we need 
to add to (}{z) a quadratic Van Vleck polynomial V{z) = V2z'^ + viz + vq with the 
undetermined coefficients V2tVi, vq which modifies the above matrix as follows: 



The operator d{z) + V{z) has a linear Stieltjes polynomial S{z) = siz + sq if and 
only if the vector (si, sq) is the left kernel of the latter matrix which leads to the 
system: 

'0 = si(Li+t-2); 

= si(Li^2 + vi) + so(Lo + V2); 

= si(Li4 + Vq) + so(Lo,i + vi); 
^0 = siLi^o + so(io,o + Vq). 

Setting si — 1 and V2 = — Li as was explained earlier we get sq = ^j^^^^^^ 
from the 2-nd equation. Substituting the obtained variables in the remaining two 
equations we get the system of two equations: 



which determines three (not necessarily distinct) pairs (ui, vq) which together with 
V2 — — ILi given us three required (not necessarily distinct) quadratic Van Vleck 
polynomials whose Stieltjes polynomials are of degree exactly 1. 

Now we finally describe the notion of natural multiplicity of a given Van Vleck 
polynomial of an operator Vi{z) used in the introduction. Let V{z) = VrZ^ + 
Vr-iz''^^ + ... + uq be some fixed Van Vleck polynomial of the Heine-Stieltjes prob- 
lem (jl.Sp . i.e. there exists a (not necessarily unique) polynomial solution S{z) of 
the equation of ()1.5p with the chosen V{z) — V{z). (Below we use notation from 
the proof of Proposition [T] ) 

Definition 2. Given a positive integer n let us define the n-th multiplicity '^n{V) 
of V{z) as the usual local algebraic multiplicity of the intersection of the (r + 1)- 
dimensional matrix pencil An{v-n Vr-i, vq) = ^o(z).„ + Vrlo + Vr-ih -I- ... -I- vqIt 
consisting of 'triangular' (ri-|- 1) x (n-|-r-|- l)-matrices with the set TM^ C TM{n + 
1, n + r + 1) of positive corank matrices at the matrix Agj^) „ -f- Vrlo + Vr-ili + 
... -|- Voir- Here (as above) A^(^z).n denotes the matrix of the action of 0(z) on the 
space Poln taken w.r.t monomial basis and Ag^^) „ -I- Vrlo + ^'r-i^i + ••■ + voir is, 
therefore, the matrix of action of the operator J)(z) -I- V{z) on Poln. In case, when 




L1-I-U2 Li^2 + vi Li^i+vo 
Lo -f U2 Lqs + vi 




(Li - ho){vo + ii,i) + {vi + Li,i)(wi + ii.2) = 0; 
{vo + io,o)(i'i + Li.,2) + (Li - Lo)ii,o = 0. 
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-^ii{z),n + 'Vrlo + Vr~ili + ... + uo-^r does not bclong to TAI^ C TM{n-\- 1, n + r + 1), 
i.e. the operator + V{z) does not annihilate any polynomial of degree at most 
n we set tJn(V^) — 0. 

Remark 13. The natural multiplicity of Van Vleck polynomials in Theorems |¥] and 
[5]while counting those with Stieltjes polynomials of degree at most n is exactly the 
n-th multiplicity from Definition [2] 

Obviously, for any given Van Vleck polynomial V{z) the sequence {tln(l^)}, n = 
0, 1, ... is a non-decreasing sequence of non-negative integers. Moreover the following 
stabilization result holds. 

Lemma 5. For any non- degenerate operator d{z) the sequence {in(V)} of multi- 
plicities of any its Van Vleck polynomial V{z) stabilizes, i.e there exists riy such 
that for all n > Uy one has jin(V^) = 'in^iV). 

Proof. Indeed, as was mentioned in e.g. the proof of Proposition [T] the leading 
coefficient Vr of V{z) must necessarily coincide with — L,„ for some non-negative m. 
The sequence { |Lj | } is strictly increasing starting from some jo , see (|1.6I) . Moreover, 
by Proposition [T] if L„ 7^ = 0, 1, n — 1 then the total multiplicity of all Van 
Vleck polynomials whose leading term equals — L„ equals ("^'^) . Therefore, if we 
take the index value jo such |Lj| > |L„j| for all j > jo then the multiplicities tJj(l^) 
can not change for j > jo since the total multiplicity increase is obtained on Van 
Vleck polynomials with a different leading coefficient when j grows. □ 

3. Proof of generalized Polya's theorems 

Let us now prove Theorem [8] following straightforwardly the recipe of [31] which 
in its turn is closely related to the proof of the classical Gauss-Lukas theorem. 

Proof. Let (zi, . . . , z„) denote the set of all roots of a Stieltjes polynomial S{z) of 
some degree n satisfying the equation (jl.Sp with a^'s being complex and /3j's being 
positive. Then for each one has 

This equation has the form 

E^+E^ = o^ (3.1) 

s=i j=i '^i 

where {^i,...,£,p) is the set of all roots of S'('''^^^(z) with p = n ~ k + 1 and 
(mi, ...,mp) is the set of multiplicities of the roots (^1, . . . ,£,p). Notice that by the 
standard Gauss-Lukas theorem all {£,1 , ■ . ■ , £,p) lie in the convex hull of the set of 
roots (21, . . . , Zn). Assume now that the convex hull of (zi, . . . , z„) is not contained 
in the convex hull of (ai, ...,«;). Then there exists some root Zj and an affinc line 
L C C separating z^ from the rest of z^'s together with all a^'s and ^m's. But then 
the equation (j3.ip can not hold since all the vectors z^ — and Zi — aj lie in the 
same half-plane. □ 

Remark 14. The above argument works for the roots of Van Vleck polynomials 
V(z) as well and extends to the case Pj > 0. 

To settle a much more delicate Theorem[S]we will prove a number of localization 
results having an independent interest. 
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3.1. Root localization for Van Vleck and Stieltjes polynomials. 



Definition 3. Given a finite (complex- valued) measure /i supported on C we call by 
its total mass the integral d^{(). The Cauchy transform C^{z) of /i is standardly 
defined as 

- I (3.2) 
Obviously, C^ (z) is analytic outside the support of /i and has a number of im- 



portant properties, e.g. that /i 



dz 



understood in the distributional sense. 



Detailed information about Cauchy transforms can be found in [16j . 

Definition 4. Given a (monic) polynomial P{z) of some degree m we associate 
with P{z) its root-counting measure lJ-p{z) = m^j^(^ ~ ^i) "^here {zi,...,Zm} 
stands for the set of all roots of P{z) with repetitions and 5{z — Zj) is the usual 
Dirac delta- function supported at Zj. 

Directly from the definition of fJ-p{z) one has that for any given polynomial P{z) 
of degree m its Cauchy transform is given by C^j,(z) = ^p^^-^ ■ 

We start with a rather simple estimate of the absolute value of the Cauchy 
transform of a probability measure which will help us to prove Proposition HI comp. 
Lemma 2 in jSj. 

Lemma 6. Let ^ be a probability measure supported in a disk Dq of radius Rq and 
centered at zq. Then for any z outside Dq one has the following estimate of the 
absolute value of its Cauchy transform C f^{z) : 

^ > \C,{z)\ > ^ 



|z - zo - i?o 



2|z - zo| ' 



(3.3) 



Proof. The l.h.s. of the above inequality is quite obvious. By (13. 2p one has that 
|C^(z)| will be maximal if one places the whole unit mass of /i at the point which 
has the least distance to z in the admissible support. In our case such a point p is 
the intersection of the boundary circle of Dq with the segment (z, zg). Its distance 
to z equals |z — zq — -Ro| which gives the required inequality. To settle the r.h.s. let 
us assume for simplicity that zq — 0. Translation invariance of our considerations 
is obvious. Let us use (j3.2|) and change the integration variable as follows: 

111 11 



C z l-Q/z 



1 



where 9 



_ C 



Since z lies outside Do and C, lies inside Dq one has |6'| < 1 which 



implies for w 



\w-l\ = 



that one has Re(i(;) > i. Indeed, 

\a\ 

= \0\\w\ < \w\ <^\w-l\< \w\ Re{w) > 



1 



1-i 



Therefore, 



dKO 



dKO 



wd^i{C) 



> 



Re{w)dn{C) 



> 



1 

□ 



Using Lemma [5] we now settle Proposition 2] (see Introduction). 



Proof. Take a pair {V{z), S{z)) where V{z) is some Van Vleck polynomial and S{z) 
is its corresponding Stieltjes polynomial of degree n. Let ^ be the root of either 
V{z) or S{z) which has the maximal modulus among all roots of the chosen V{z) 
and S{z). We want to show that there exists a radius R > such that |^| < R for 
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any ^ as above and as soon as n is large enough. Substituting V{z), S{z), ^ in (jl.Sp 
and using (|1.4p we get the relation: 



dividing which by its first term we obtain: 

Notice that the rational function :— , ^ -sgf-f/ x 

the polynomial S'(*^(z). Easy arithmetic shows that 



(3.4) 



is the Cauchy transform of 



(n-A:)! 



in-k + ~ z) n b,iz) S(>^)(z) _ ,)! n^^- &,(^) 



Notice additionally, that by the usual Gauss-Lucas theorem all roots of any S^^'^ (z) 
lie within the convex hull of the set of roots of S{z). In particular, all these roots 
lie within the disk of radius |^|. Therefore, using Lemma [6] we get 



(3.5) 



Notice that since Qk{z) is a monic polynomial of degree fc + r (recall that r is the 
Fuchs index of the operator f (z)) then one can choose a radius R such that for any 

z with \z\ > R one has \Qk{z)\ > — . Now since for any i = 1, k — 1 one has 
degQi{z) < i + r we can choose a positive constant K such that |Qi(z)| < i^lz]'^'" 
for all i = 1, ...,k — 1 and \z\ > R. We want to show that can not be too large for 
a sufficiently large n. Using our previous assumptions and assuming additionally 
that 1^1 > i? we get 



\Q,iO\ (n-fc)! ,fe_, K ■ 2^-^+1 



i)...{n -k + l)' 



Now we can finally choose A^o large enough such that for all n > A'o, all i 
and any |^| > i? one has that 



,fc-l 



But then obviously the relation 



< 



K ■ 2 



k~i+l 



< 



1 



{n - i)...{n - k + 1) k - 1 

can not hold for all n > Nq and any |^| > iZ 



E 



fc-i 



k-l 

<E 



1 



< L 



□ 



Now we will strengthen the arguments in the proof of Proposition U] in order to 
settle Theorem El Denote by TZq,, the maximal distance between the origin and 
ConvQi^. The following statement holds. 

Lemma 7. For any non-degenerate higher Lame operator d{z) and a given number 
S > there exists a positive integer Ng such that the roots of all Van Vleck poly- 
nomials V{z) possessing a Stieltjes polynomial S{z) of degree > Ng as well as the 
roots of these Stieltjes polynomials lie in the disk \z\ < TZq^ + S. 
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Proof. Notice that once S is fixed the quotient jg^^^'^^^i is bounded from above for 
each i = 1, A: — 1 if we assume that \z\ > TZq^ + S. Indeed, all the roots of 
Qk{z) lie within the disk of radius TZq^ centered at the origin and each Qi{z) has 
a smaller degree than Qk{z). Consider now again the estimate (j3.5p . Since we now 
know that ^ lies in some bounded domain for all possible polynomials V{z) and 
S{z) of sufficiently high degree and that the quotient is bounded from above 

outside the disk of radius TZq^ + (5 we get that the right-hand side of p.Sp goes to 
when n — > oo under the assumption that ^ stays outside the latter disk. Looking 
again at (j3.4p we see that by the latter argument it can not hold for |^| > TZq^. + 5 
when n — > oo. This contradiction proves the lemma. □ 

To finish the proof of Theorem [9] notice that the choice of the origin is in our 
hands, i.e we can make an arbitrary affine shift of the independent variable z and use 
the same arguments. Since the convex hull Convq^ is the intersection of all disks 
centered at different points and containing Convq^ we can for any chosen e > find 
the intersection K, of finitely many disks in C such that /C contains Convq^. but is 
contained in Couvq^ . (One can choose one such disk for each edge of the boundary 
of ConvQ^. putting its center sufficiently far away on the line perpendicular to the 
edge and passing through its middle point.) Then since /C is the intersection of 
finitely many disks we can applying Lemma[7]find such that all roots of all V{z) 
and S(z) for all n > iVf lie in /C. □ 

4. 'On the existence and number of Lame functions of higher 



4.1. Comments on Heine's result and history around it. As Heine himself 
mentions in [18 the requirement of algebraic independence of the coefficients of 
(52 (-z) and Qi{z) being sufficient is too strong and restrictive for his purposes but 
he fails to give any other explicit condition guaranteeing the same result, see The- 
orem [TJ Heine's original motivation for the consideration (II. 2p comes from the 
classical Lame equation (II. ip in which case Qi{z) and Q2{z) are very much alge- 
braically dependent, namely, Qi(z) = Q^iz)!^. Also in order to prove that the 



Qi{z) Heine uses an inductive argument where he forces Qi{z^ and Q\{z) to become 
algebraically dependent in a special way. Theorem[2]is another clear indication that 
the algebraic independence is apparently an inappropriate condition for the goal. 

Interpretation of Heine's text written in a rather cumbersome 19-th century 
German and the exact statements it contains seems to create difficulties for math- 
ematicians starting from 1870's and up to now, see e.g. [37j. The main classical 
sources, namely, jH] from 1885, [3T] from 1912 and [12] from 1939 are not too clear 
about what it is that Heine actually proved and under what assumptions on Q2(^z) 
and one can guarantee that for a given positive n the number of possible 

polynomial pairs (y, S) such that the corresponding S has degree exactly n is finite 
and bounded by ("+^"^). Being aware of the existence of a gap in his proof Heine 
seems to covers by a reference to a letter of his friend Leopold Kronecker who has 
(under unknown conditions) shown that for a given degree n the eliminant of the 
system of algebraic equations defining the coefficients of the polynomial Viz) does 
not vanish identically. This statement is equivalent to the finiteness of the number 
of these polynomials. Heine mentions also a short note of Kronecker's on this topic 
presented in the January issue of Monatsbericht der Berliner Akademie from 1864. 
Unfortunately the track ends here. All one can find in this issue is the phrase that 
at the meeting on the section of physics and mathematics of the Prussian Academy 
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upper bound 
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of Science held on November 14, 1864 "Herr Kronecker las iiber die verschieden 
Faktoren des Discriminante von Eliminantions-Gleichungen ", i.e. "Herr Kronecker 
gave a lecture on different factors of the discriminant of the elimination equation." 
In attempt to overbridge this gap we undertook the task of translating and (even 
more so) decoding Heine's arguments. In what follows we give a frase-by-frase 
translation of portions of §135 and §136 from Heine's book [TH] which are relevant 
for our consideration. Some of the phrases were difficult to understand literally 
and we provided our interpretation and comments of the content placed within the 
slash signs. We allowed ourselves to correct several obvious misprints without a 
special mention, tried to keep (to certain extent) the language flavor and preserved 
the enumeration of the formulas of the original text. 

4.2. Translation. §135 ... 

Next we ask which conditions the polynomials xi^) 9{x) must satisfy in 
order for the differential equation 

d'^W dW 
^{x)^+x{x)-g+e{x)W = (88) 

to have a solution which is a polynomial of degree n in a; assuming that is of 
degree p + 1 and x ^^id 9 are of degrees at most p resp. p — 1. /In fact, x(a;) is 
fixed and 6{x) is a variable polynomial. It is not clear why Heine talks about both 
x{x) and 9{x) here./ 

This is always the case when it is the question about, as in the case of Lame 
functions, a differential equation whose general solution does not contain any higher 
transcendentals than a rational function of integrals of algebraic functions and 
which has a definite order at a; = oo. /This and the next two phrases explain why 
one should assume that deg(x) < deg(V')./ 

We say about a function W that it is of the order a for the finite value a; = a, if 
{x — a)°'~'^W and {x — a)°'^'^W will be respectively oo for x = a, however small e 
is taken; we assign to it the order a at infinity, if x~°'~^W and x~°'^'^W for x — oo 
becomes resp. oo. Thus, for example, logx has a definite order, namely 0. 

If y and z are two particular solutions of (88) then one has 

log(yz' - zy') = J ^dx. 

If X were not of smaller degree than ip then yz' — zy' would at a; = oo go to or 
oo as an exponential function and would therefore have no order. 

For the solution to have an order for every x where tpix) vanishes, after 
possible cancellations can in the denominator only have different factors. This 
follows from the same equality between two particular solutions which we used 
above. 

The following theorem answers the question posed at the begining: 

// the two polynomials ip{x) and x(a;) are given, the first of degree p + 1, the 
second of degree p, then for exactly y~ i) different functions 9{x), 

there exists a particular solution of (88) which is a polynomial of degree n in x. 

For p = 1 we understand the number given above which in general might be 
denoted by (n,p) as 1. It is assumed above that the coefficients of '0 and x are 
mutually independent. I call the numbers a, 6, etc. mutually independent when 
there is no algebraic equation with integer coefficients which they satisfy. /As stated 
above Heine actually proves his theorem under the assumption of the algebraic 
independence of the coefficients./ 
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It may be immediately added here that we will say about the numbers a,b, . . . 
which already satisfy one or more algebraic equations that "they are further spe- 
cialized" when in addition to these equations they satisfy one or more additional 
equations - which of course are not allowed to contradict the earlier ones. 

The just mentioned assumption for the validity of the theorem demands more 
than is necessary. One could say with the same right that a polynomial of degree n in 
X with mutually independent coefficients vanishes for n different values of x while for 
this it would suffice that the coefficients did not satisfy a particular equation with 
integer coefficients, namely, the well-known one which describes the coincidence 
of roots. Also for the validity of our theorem it suffices that the coefficients do 
not satisfy certain finite number of algebraic equations, which equations in every 
case can be found but not in a comprehensible way. /As we mentioned in the 
introduction Heine realizes that for a given fixed n the set of all pairs of polynomials 
{ip{x),x{x)) of degrees at most p + 1 and p respectively for which there are less 
than (n,p) functions 0{x) solving the problem under consideration is an algebraic 
hypersurface with integer coefficients. It is by no means clear to us how Heine 
could possibly conclude that all the coefficients of the discriminantal equation are 
integers in the basis of the coefficients of ip{x) and xi^)- The development of the 
corresponding theory can be traced back to Cay ley, see appendix in 17 \, but no 
general results were obtained until much later./ 

§136. To obtain the proof of the theorem one substitutes in (88) polynomials of 
degree n and p — 1 for W and 9, namely 

W^ = x"+gia;"-i + ...+, 

9 = koxP-^ + kixP-^ + k2xP-^ + ... 

It is clear that the necessary and sufficient condition that W satisfies the equation 
(88) is that certain n + p equations are satisfies which are linear in both the ^^'s 
and fcj's and in the coefficients of V' and x- To show the structure of these, without 
having to work with too clumsy formulas I present them for the case p = 3. 

Let the given functions be 

tpix) = x{coX^ + CiX^ + C2X + C3), 

xix) = box^ + hix^ + 622; + &3 

and the sought functions 

W{x) = 50a;" + gia;""^ + g2x"-=^ + . . . + .g„, 

9[x) — kox'^ + kix + k2- 

/Notice that Heine apparently realizes that, in general, it might be impossible to 
find W{x) as a polynomial of degree exactly n and introduces even the leading 
coefficient as a new variable without special explanations./ 

For W to satisfy the differential equation (88) the coefficients giS and fcj's must 
satisfy the system of equations: 
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=5o[fco + nbf) + n{n - l)co]; 

=5i[fco + (n - l)feo + {n- l)(n - 2)co] +go[fci + ^i&i + n{n - l)ci]; 

=g2[ko + (n - 2)6o + (n - 2)(n - 3)co] + ffi[/ci + (n - + (« - l)(n - 2)ci] + 

+ 50 [^2 + nb2 + n{n - l)c2]; 
=gs[ko + {n- 3)6o + (n - 3)(n - 4)co] + 5(2[A;i + (n - 2)6i + (n - 2)(n - 3)ci]+ 

+ fliffe + - l)b2 + (n - l)(n - 2)C2] + go Ha + n{n - 1)03]: 
=5f4[A;o + (n - 4)60 + - 4)(n - 5)co] + 33(^1 + - 3)6i + (n - 3)(n - 4)ci] + 

+ 92[k2 + (n - 2)62 + (n - 2)(n - 3)c2] + <?i[(n - 1)63 + (n - l)(n - 2)03]; 

In this way the equations wiU continue to be formed so that the next one wih 
give the relation between the four gfj's with the indices 5, 4, 3, 2. The final equations 
will be 

=gn-i[ko + ho] +gn-2[ki + 26i + 2 ■ Ici] + ff„-3[fc2 + 36i + 3 • 2c2] +ff„-4[463 + 4 
0=5„[fco] +5f„_i[fci + l-6i]+5f„_2[fc2 + 262 + 2-lc2] +fif„_3[363 + 3 

0= gn[ko]+9n-l[kl + l-hi] +fi(„_2[263 + 2 

= 9n[ko] +gn- 

From the first equation the coefficient fco is completely determined in terms of the 
given coefficients 60 and cq of ip and x; the next n equations give all gi's expressed 
in terms of the same known coefficients 6i's and c^-'s and the (p— 1) (in our example 
2) unknowns fci, fc2, . . .. The values of g^'s that are obtained from the second to 
the (n + l)st equations when substituted in the last (p — 1) equations, will then 
give {p — 1) equations of higher degrees between the unknowns fci, ^2, • • • , kp-i and 
the known coefficients of ijj and x^ which only appear rationally in these equations. 
/Heine apparently means that g^'s will be given by rational functions of c;'s and 
6to's and substituting these one gets a system of rational equations defining fc^'s./ 

Once the fcj's have been determined from these p — 1 equations, the substitution 
of the found values in the 2nd to the (n + l)st equations will give the gi's. Two 
systems of related fcj's are called different, i.e. the two systems fci, . . . , fcp_i and 
k[, . . . , kp_i are called different when they are not equal. One realizes with the full 
confidence from the form of the 2nd to the (n+l)-st equation that every system fc^'s 
corresponds to a system of gt's and different systems of kj's correspond to different 
system of g^'s. /The latter statement of Heine is false as is. It requires that the 
diagonal entries in the uppertriangular system are non-vanishing (see our nonres- 
onance condition in the introduction). But it is certainly true if the coefficients of 
ip{x) and x(a;) are algebraically independent./ 

One obtains thus that 
There exist as many different equations (88) and therefore as many different poly- 
nomials W 0/ degree n as there are different systems of kj 's. 

Next it is realized that 
The degree of the elimination equation is at most {n,p), that is there can only be 
at m,ost [n.p] different systems of kj 's. 

If one throws a glance at the (n + p) equations, which, with the exception of 
the first one for ko, one can find above for the special case that p = 3, one will 
perhaps not realize the truth of this statement immediately, and instead believe 
that the degree of the elimination equation is larger. / One is supposed to disregard 
the 1st equation in the system above and study the remaining (n + p) equations./ 
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But if one instead of fc2, fcs, . . . , substitutes X2,x^, . . . where the lower numbers are 
indices and the upper numbers are exponents, and for symmetry we use Xi ki, 
one reahzes immediately that gi, §2, . ■ . , g-n are polynomials in the XiS of degree 
1, 2, . . . , n resp. so that after the substitution the (p — 1) last equations will have 
degrees n + 1, n + 2, . . . , n + p + 1 in x^'s. That means that the degree of the 
elimination equation will at most grow to (n + l)(n + 2) . . . (n + p + 1). If one 
takes into account that every value of ki,k2, ■ ■ ■ corresponds to one value of X\, 
two values of X2, three values of x^ then the above assertion is proved. /This is an 
excellent passage! What Heine does is called in the modern language of algebraic 
geometry the weighted Bezout theorem, see e.g. [13]. The author was unable to 
find a reliable proof of this result in the literature prior to 1970./ 
Under the assumptions that the elimination equation is not identically zero it will 
indeed have the above mentioned degree and give {n,p) different systems ofk's. 
/Crucial claim but not completely proved below./ 

There exist as I will show below indeed (n, p) distinct systems of the coefficients 
if the coefficients of tp and x are specialized in a certain way. /Heine will show 
by induction that for a special choice of "0 and x oii6 can obtain exactly (n,p) 
distinct solutions. But instead of algebraically independent coefficients of and x 
he needs to make them dependent to get an example of {n,p) distinct solutions. 
This is correct as soon as one knows that even for the specialized situation the 
total number of solutions is finite. This finiteness probably follows from his specific 
choice of specialization, see below./ 

That this elimination resultant does not vanish identically follows from the next 
observation which I take from a letter of my friend Kronecker. /One needs to check 
also that the system of equations has (under very unclear non-degeneracy assump- 
tions on the coefficients of ip and x) only isolated solutions. This is equivalent to 
the non- vanishing of the eliminant./ 

If in the mentioned final equation which will determine the functions 0{x) and 
W{x) all coefficients vanish then according to the general principle of elimination 
will at least one of the roots of W{x) = be unrestricted. /This is an interesting 
although a rather obvious observation. Notice that the elimination theory hardly 
at all existed in 1870's./ 

If one assigns to this root all values for which 4>{x) vanishes then one gets through 
this (procedure) certain restrictions on the function x(^) but the latter function 
does not satisfy them even after the mentioned specialization. 

Herr Kronecker added in the mentioned message that these restrictions are ac- 
tually satisfied and one of the roots of W{x) =0 remains undetermined if both 
ijj{x) and x{^) have the properties that for the known function 9(x) both solutions 
of (88) are polynomials in x. 

Concerning the number of systems with specialized ip and x I set such relations 
between the coefficients which give that ^ has a factor {x — a) twice and x has 
it once. /Then the algebraic independence is lost here since ip has a double root 
and therefore lies on a discriminantal surface./ Then all W satisfying (88) have the 

^ In Monatsbericht der Berliner Akadcmie from January 1864 added (ackomplished) Herr 
Kronecker my message with the following 

Introducing the roots of W{x) = as variables in the equations 

V'(xfc)VK"(xfe) + xi^k)W'(xk) = 0, fc = 1, 2, . . . ,n 

which define them and substituting the coefficients of W and W" through the symmetric functions 
o{ xi, . . . ,Xn one sees directly that one of the unknowns x remains arbitrary if the elimination 
equation vanishes. Through a simple transformation of this system of equations one can determine 
the degree of the final equation and at the same time prove that certain coefficients are different 
from as long as there are no special conditions on the functions ip and x- /A rathe unclear proof 
of the nonvanishing of the eliminant./ 
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form 

U{n); {x - a)U{n - 1); . . . ; (x - a)"C/(0), 
where U are as in §123 polynomials coprime with {x — a) are their degrees are given 
in parenthesis to the left of the letter U. Under the substitution of this expression 
in (88) one gets for every U an equation like (88) in which instead of ip and x 
appear polynomials with unrestricted coefficients which arc not of degree p+l and 
p but instead of p and p — 1. /Apparently the logics here is as follows. We can use 
our result to prove that the number {n,p) of simple solutions is obtained if we can 
show that for the above specialization the total number of solutions is finite./ 

If one assumes that the general statement which we are proving is settled if ip is 
the product of p linear factors (and for the product of 2 linear factors this is easy 
to show) then one gets the situation when consists of p + 1 factors of which two 
are coinciding, alltogether 

(n,p - 1) + (n - l,p - 1) + (n - 2,p - 1) + . . . + (0,p - 1) 

which after summation gives {n,p) different VF, i.e. (n,p) different 9 just as many 
as different systems of fc's. /This accomplished the induction step. What one 
misses is mentioning that the total number of solutions for Heine's specialization is 
finite. But this follows from his representation of all solutions as U{n); {x — a)U{n — 
1); . . . ; {x — a)"U (0) ,. In each of these cases we already know that the number of 
solutions is finite and all of them are simple. To be completely rigorous one should 
use a double induction on p and n as we did in § [5] Additional simplification of 
the order 2 case compared to the general order k case considered above comes from 
the fact that during the above specialization one can assume that the polynomials 
ip and X such that ip = ~ a)'0 and x = — a)x have algebraically independent 
coefficients which does not work for higher order case./ 

5. Final Remarks 
Let us formulate a number of relevant questions and conjectures. 

Problem 1. Is it possible to describe when a linear ordinary differential equation 
with polynomial coefficients admits at least 2 polynomial solutions? 

The prototype result of Varchenko-Scherbak gives a satisfactory answer for equa- 
tions of the second order. The answer to the latter question allows to detect the 
appearance of multi-dimensional families of Stieltjes polynomials. 

Problem 2. Under the nonresonance assumption (II. 7p is it possible to obtain 
explicitly the discriminantal surface which shows when a Van Vleck polynomial 
attains a non-trivial multiplicity. 

This question addresses the problem of explicit determination of the discrimi- 
nantal surface mentioned in Heine's proof. Some discussion of this problem can be 
found in [39] . 

Problem 3. Explain how the number of Van Vlecks polynomials having Stieltjes 
polynomials for a certain given degree n can drop below ("^'^) ? 

The next question addresses the issue of location of the roots of Van Vleck and 
Stieltjes polynomials. 

Problem 4. Under what assumptions on t)(z) the roots of any its Van Vleck and 
Stieltjes polynomials lie in the convex hull of its leading coefficient Qk{z)l 
The basic examples are provided by Stieltjes's and Polya's theorems. 
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Finally, 

Problem 5. Is it possible to extend the results of this paper to the case of degen- 
erate higher Lame operators? 

T. Bergkvist [5] has obtained a number of interesting results and conjectures 
in the case of degenerate exactly solvable operators. Motivated by her results we 
formulate the following conjecture. 

Conjecture 1. For any degenerate Lame operator and any positive integer Nq the 
union of all the roots to polynomials V and S taken over degS" > Nq is always 
unbounded. Therefore, this property is a key distinction between non- degenerate 
and degenerate Lame operators. 
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